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Abstract 

The  equations  of  motion  of  a  ship  constrained  to  oscillate  in 
heave,  pitch,  and  roll,  while  moving  ahead  with  a  uniform  velocity 
in  oblique  waves  are  studied.  It  is  shown  that  because  of  the  non¬ 
linear  static  coupling  which  is  known  to  exist  between  theqe  degrees 
of  freedom,  the  rolling  response  can  became  large  when  the  period  of 
wave  encounter  is  in  the  neighborhood  of  one  half  the  natural  rolling 
period.  In  this  case  the  response  has  a  period  which  is  equal  to 
twice  the  period  of  the  excitation,  and  which,  as  a  consequence,  is 
nearly  equal  to  the  natural  period. 


ii 


U&L  of  ^TmbolS 

x,  y»  *  *  displacements  of  fhip  along  x,  y,  z~axas  (surge, 
ayay  end  he eve) . 

^  0^  y!  m  rotations  about  x,  y,  z-exis  (roll,  pltc}i  and 
y aw) . 

u  r  component  of  ship  velocity  along  x-axes. 

j  T  *  mass  moment*  of  Inertia  about  x,  y-exes. 

*  -1 

^  M  r  components  of  moment  along  x,  y-axes. 
f  component  of  force  along  z-axis. 

2A  «  wave  height. 

L  -  ship  length. 

A  ■  wave  length. 

60  *  frequency  of  wave  encounter. 

0{  -  heading  of  ship  with  respect  to  direction  of 
wave  travel. 

^  -  fraction  of  critical  damping. 


Dots  denote  time  differentiation. 
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Introduction 


In  the  study  of  ship  motions,  ss  in  ths  study  of  ths  motion  of 
spy  rigid  body,  two  approaches  are  generally  used  to  refine  the 
llnf it  theory:  One  may  Introduce  refinements  In  the  terms  already 
''.-appearing,  in  the  linearised  equations  of  motion  (e.g.  added  mass 
and  damping  terms)  or  one  may  Introduce  higher  order  terms  into  the 
equations.  It  is  the  latter  approach  which  is  adopted  in  this  paper. 

,  *•  '•**..  '  **  i 

r’wj&k,  nonlinear  terms  retained  in  the  equations  of  motion  are , 
part,  the  same  second  order  static  coupling  terms 
Paul ling  and  Ros.enberg  [l]1,  Kinney  [R,5]>  ’end  Hsu  [3,4], 


The;., difference  between  the  investigation  presented  here  and  the  ones 

jiist  cited  is  that  here  we  consider  the  ship  free  to  oscillate  simul* 

*  ,*  -  **'*»• 

•  .  ■  t  . 

taneously  in  the  degrees  of  freedom  of  heave,  pitch,  and  roll  while 
.mdving  ahead  with  uniform  velocity  in  oblique  waves,  and  Include  the 
effect 'b£  second  order  damping  in  roll.  It  Is  shown  that  when  the 
frequency  of  wave  encounter  is  in  the  neighborhood  of  twice  the 
natural  rolling  frequency,  the  roll  amplitude  may  tend  to  become 
large,  and  the  rolling  response  may  exhibit  a  period  equal  to  twice 
the  period  of  the  excitation,  which,  as  a  consequence,  is  nearly 
equal  to  the  natural  period. 


Hiasbers  in  brackets  designate  References  at  end  of  paper. 
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1.  Equations  of  Motion 

Here,  we  shall  derive  the  equations  o£  motion  of  a  rigid  ship 
executing  motions  about  Its  equilibrium  position.  The  purpose  of 
this  discussion  Is  two-fold: 

(a)  To  state  clearly,  and  to  justify,  all  assumptions  considered 
necessary  to  render  the  problem  tractable,  and 

(b)  To  demonstrate  In  what  way  all  earlier  analyses  [1,2, 3, 4, 5] 
constitute  special  cases  of  this  paper. 

The  coordinate  axes  are  the  xyz-trlad  whose  origin  coincides 
with  the  mass  center  of  the  ship.  The  x-axis  lies  In  the  fore-and- 
aft  plane  of  symmetry  of  the  ship;  It  Is  parallel  to  the  design 
water  plane  of  the  ship  and  Is  positive  forward.  The  y-axls  Is  posi¬ 
tive  to  starboard,  its  eleveation  being  such  that  the  xy-plane  is 
parallel  to  the  design  water  plane.  The  z-axis  Is  orthogonal  to  the 
x  and  y-axes,  and  Is  positive  downward.  The  six  degrees  of  freedom 
of  the  ship  consist  of  translations  parallel  to  the  xyz-axes,  and 
rotations  about  them.  The  translations  will  be  denoted  by  x,y,z 
respectively,  and  the  rotations  by:  <fi  about  the  x-axls,  Q  about 
the  y-axls,  and  about  the  z-axls.  The  translations  are  positive 
in  the  positive  directions  of  the  axes,  and  the  rotations  are  posi¬ 
tive  in  accordance  with  the  right-hand  rule.  All  coordinates  are 
sero  in  the  equilibrium  position. 


The  ship  is  rigidly  constrained  in  sway  and  yaw,  l.e.  js^sO  , 
end  the  surging  motion  Is  prescribed  to  be  uniform. 


Justification: 

The  above  assumption  is  made  because  it  Is  the  explicit  purpose 
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of  this  analysis  to  examine  only  motion  In  the  remaining 
degrees  of  freedom. 

Aasvaption  B: 

The  xyz-axes  are  principal  axes  of  the  ship. 

Justification: 

Ships  to  be  considered  here  possess  athwertshlp  symmetry  and 
nearly  fore-and-aft  symmetry.  The  deviation  from  the  latter 
is  so  small  that  this  assumption  is  considered  justified. 

Under  Assumptions  A  and  B,  the  equations  of  motion  are  those 
in  heave,  pitch,  and  roll;  they  are 

- 

_ I u-6  *  M - (l.i) 

- I^Ls.  R  - 

where  m  is  the  mass  of  the  ship,  Iy  and  lx  are  the  mass  moments  of 
inertia  about  the  y  and  x-axes  respectively,  and  ?  )  ft  ,  and  K 
are  the  external  forces  or  moments  in  heave,  pitch,  and  roll 
respectively. 

Assumption  C: 

Z  ,  M  ,  and  k  are,  in  the  neighborhood  of  the  equilibrium 
position,  analytic  functions  of  the  displacesients  96  >  9  , 
the  velocities  ^  ,  <^>  ,  6  ,  and  the  accelerations  ^  ,  9  . 

Justification: 

If  Assumption  C  were  violated,  linear  analysis  about  the  equi¬ 
libria  would  be  lnadkalssable . 

Thus,  the  right-hand  sides  of  (1.1)  may  be  represented  by  a  general 
expression 
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where  *2 \ G2,9  ft  )Qj«  1C  i  and  since  the  coordinates  vanish  in  the 

equilibrium  position, 

Q.(o,o,o,o,o,olo,01o')  =  0 

Moreover,  in  view  of  Assumption  C,  every  Qj  may  be  expanded  in  a 
Taylor's  series  about  the  equilibrium  position. 

Assumption  D: 

The  Taylor's  series  expansion  of  the  Q*  will  Include  only 
terms  up  to  and  including  second  order. 

Justification: 

It  is  the  explicit  purpose  of  this  analysis  to  examine  certain 
second  order  effects. 

Writing  temporarily  ^  ^  <p  »  the  terms  on  the 

right-hand  side  of  (1.1)  became 

x  y  t  4  ^  \ 


j-i 

+z  e*.  £■  ■ 4 


(1.2) 


+  ^  \  i  'i  |  ii.TL  \ 


where  the  partial  derivatives  are  evaluated  at  the  equilibrium 
position. 
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As  sumption  E: 

All  second  order  terms  In  (1.2)  Involving  velocities  end 
ecceleretlons  (except  for  the  term 
discarded. 

Justification: 

This  assumption  is  dictated  largely  by  the  lack  of  knowledge 
of  the  magnitude  of  the  ignored  derivatives.  The  only  term 
retained  here  is  the  second  order  damping  effect  in  roll;  its 
physical  meaning  is  clear,  its  magnitude  can  easily  be  deter¬ 
mined,  and  its  effect  is  considered  significant.  There  are, 
however,  other  rational  reasons  for  this  assiaaptlon.  Since  the 
present  analysis  is  a  refinement  on  linear  analysis,  one  may 
safely  assume  that  the  primary  effects  of  velocity  and  accel¬ 
eration  on  the  motion  are  Included  by  retaining  the  first  order 
terms  in  these  quantities.  Moreoever,  this  analysis  is  con¬ 
cerned  with  the  effect  of  static  coupling  on  the  stability  of 
the  motion,  and  the  discarded  terms  do  not  represent  static 
coupling. 

Under  all  assumptions  listed  here,  the  equations  of  motion 

became 

i  -  6  6  +l\i**s6 

♦***+**$♦*  y>  a.3) 

$  =  M(  e  +  +  M**  * 

m*  *i  m„  *  r  *  2  m,/1 + Mi,»  \6 

+ 8  +  ]s  • 


)  will  be 


(1.4) 
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0  *  Kp <P  +  ^ +  +  ^  ^  0 

+  <g  6  +  K..  j  +  £  K^'  ♦  i  K,/  <- 1  ^  ^ 

+w*  v>fV*  ■ 

where  Che  subscripts  denote  pertlel  differentietion,  where 

'etc>  V-A — t+c' 


(1.5) 


£  _  - ?Jl_ 

B  -  2  •  • 


K  ^ 


(1.6) 


^  etc. 


end  where  in  (1.5)  the  term  which  would  normally  appear  la 

the  Taylor's  expansion  has  been  replaced  by  The 

reason  for  this  is  that  physfoally  the  term  represents  a  damping 
maaMnt  and  therefore  should  change  sign  with  the  velocity. 

The  equations  (1.3),  (1.4),  and  (1.5)  simplify  considerably 
when  the  athwart ship  symmetry  of  the  ship  is  utilized.  In  particular 

W  W‘z*"° 
WVVM«*° 

K  K»  *  <  V *  K»  ’  K  \  ’ ' **'  ■  K*p  Ke«“  K\\  ke  *  0 

so  that  the  equations  of  motion  become 

a.?) 

ji=  M4e*M.s  +  MvT5+Mi\  +  M^'^ 


(1.8) 
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t-w+w*  +K+>*6- 


(1.9) 


Next  we  propose 
Assumption  F: 


Z  =  M  =  0 
ii 


Justification: 

It  can  readily  be  shown  numerically  that  the  linear  terms  ^ 
and  Mg  6  account  adequately  for  the  heaving  force  due  to 
heaving,  and  the  rolling  moment  due  to  roll.  Moreoever,  the 
second  order  terms  in  this  assumption  do  not  represent  static 
coupling  between  degrees  of  freedom. 

Assumption  G: 

The  effect  of  a  seaway  on  the  motion  will  be  taken  into  account 
by  adding  terms  2(t)  ,  ,  and  K.(i')  ,  respectively  to 

the  equations  of  motion. 

Jv*mU9itlqn: 

The  waves  are  considered  as  known  functions  of  time,  and  the 
effect  of  the  ship  itself  on  the  seaway  is  considered  negligible 
Under  Assunptions  A  to  G,  the  equations  of  motion  become 


i**iV*\y*>*+*e6+*s6 

e=V+Me®+MU+MiVMn 

?  -  *  W  +  K*fb  +  <W®+  <(i)  • 


(1.10) 


(1.11) 


n 


(1.12) 
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Not  all  of  the  partial  derivatives  in  the  equations  of  motion  are 
independent.  In  fact  they  satisfy  the  following  relationships  [1,2] 

^ e  =  ^  =  1  *ee = 

a  7  Ta  lx  U-13) 

Equations  (1.10),  (1.11),  and  (1.12),  together  with  (1.13)  are  those 
on  which  all  future  development  is  based. 

Below,  we  examine  the  relation  between  these  equations  of 
motion  and  those  of  earlier  analyses.  All  contained 
Assumption  A* : 


2se=  4*  “  M*j 


M 


Under  this  assumption  the  stability  analysis  of  the  motion  is 
reduced  to  the  discussion  of  Mathieu  equations  which  can  be  decoupled. 

The  equations  of  Hsu  [3]  emerge  by  considering  the  ship  as 
symmetric  fore-and-aft,  or 
AMyEtjLqnJ1.: 


^B  ~^9  ^9 


i  a 


My 

and  that  the  waves  have  sinusoidal  shape  and  travel  with  uniform 
velocity,  or 

Assumption  C  * :  ;?  ("£)  =  Ze  0  COS  oJ't 

M(t)*  Heo  sin  u)t 

where  u)  is  the  frequency  of  wave  encounter. 

In  that  case  one  finds  Hsu's  [3]  equations 


c-  \ =  ?«<>  ^ 


0  +  C Weos'*  ^ 

=0 


(1.14) 

(1.15) 

(1.16) 


where  use  has  been  made  of  the  notation 
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> 

£ 

1 

II 

& 

•« 

i 

ii 

4) 

O 

(1.17) 

> 

6  ng 

(1.18) 

Hsu's  equation 

[4]  are  obtained  by  admitting  as 

the  only  degrees 

of  freedom  those  of  heave  and  pitch,  or 

Assumption  B"  :  2^*£d  ■  2^  =  2^.  *  Mj  =  O 

and  by  assisaing  a  calm  sea,  or 
UUMBSlSB.  C"i,  2(«=  MttVO. 


In  that  case  one  finds  [4] 

If  one  also  proposes 
Assumption  JD^:  £  =  MQ 


(1.19) 


Hsu  shows  that  the  solutions  of  (1.19)  are  always  unstable. 

If  one  adults  only  the  degrees  of  freedom  of  pitch  and  roll 
and  Includes  linear  damping  In  roll,  but  no  damping  In  pitch,  one 

has 

AfS«aptV>n  B Mt:  ^  ^  =  M  (O  -  K  CO  =  0 


and  the  equations  of  motion  are  those  examined  by  Kinney  [5].  They 
are 

e  -V 


(1.20) 
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If  all  damping  is  ignored,  and  any  two  of  the  three  degrees  of 
freedom  of  heave,  pitch,  and  roll  are  a<fenltted,  one  finds  the 
equations  of  Paulling  and  Rosenberg  [1]  for 

Heave-Roll:  ..  _ 

Pitch-Roll:  .•  .  . 

Heave-Pitch:  _  _  _ 

5  «M,a 

These  last  three  pairs  of  equations  [1]  are  those  which  have 
initiated  this  study  of  the  nonlinear  static  coupling  effects  on  the 
stability  of  ship  motions. 

To  complete  this  discussion  we  shall  present  below  a  brief 
summary  of  the  results  of  the  earlier  analyses  [1,2, 3, 4, 5].  We  will 
show  that  all  the  equations  of  motion  considered  In  these  analyses 
can  be  reduced  to  the  standard  form  of  Mathieu's  equation,  the 
stability  of  which  can  be  discussed  in  a  straight  forward  manner. 

We  will  therefore  begin  by  giving  a  brief  discussion  of  Mathieu's 
equation. 

The  standard  form  of  the  Mathieu  equation  with  damping  is 


where  f<  ,  %  and  £  are  three  arbitrary  parameters . 


(1.21) 


(1.22) 


(1.23) 


(1.24) 
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The  stability  of  a  solution  y«y(x)  can  be  deduced  from  the 
Ince-Strutt  chart  which  consists  of  a  series  of  curves  which 

divide  the  $£  -plane  into  stable  and  unstable  regions.  For  values 
of  $  and  6  which  define  a  point  inside  of,  or  on  the  boundary  of 
an  unstable  region,  the  solution  is  unstable.  Figure  1  shows  an 
approximation  of  the  curves  defining  the  boundary  between  the  stable 
and  unstable  regions  for  k  ■  0,  due  to  Stoker  [6],  which  is  adequate 
for  moderate  values  of  £ 

For  lc  ykO  the  equation  of  the  boundary  curves  between  the 
stable  and  unstable  regions  near  %  m  1/4  have  been  given  in  [5] 
for  small  values  of  £  .  It  is* 

+  (1‘25 

where  0  is  Che  fraction  of  critical  damping  and  is  plotted  in 
Figure  2. 

For  a  more  detailed  discussion  of  the  solutions  to  Mathleu's 
equations  the  reader  is  referred  to  McLachlan  [7]. 

We  shall  now  examine  the  stability  of  the  solutions  of  the 
equations  of  motion  presented  above. 


2  The  standard  form  of  Mathleu's  equation  in  [5]  is  ^+2 K CS-Ze^as Z%\ 

,71*0.  instead  of  (1.24)  so  that  the  equation  of  the  boundary 
curvb  presented  there  is  in  terms  of  £.  and  a  .  The  conversion 
to  %  and  e.  is  self-evident.  0 
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FIG.  I  STOKERS  APPROXIMATION  OF  STABLE  AND 
UNSTABLE  REGIONS  FOR  THE  MATHIEU 
EQUATION 


8 


FIG.  2  THE  EFFECT  OF  DAMPING  ON  THE  &€  PLANE 
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The  equations  (1.14),  (1.15),  and  (1.16)  o£  Hsu  [3]  may  be  put 
In  standard  £orm  as  follows:  Equation  (1.14)  can  be  easily  inte¬ 
grated,  giving 


^  *■  A  cos  (tot  - 


(1.26) 


where 


V 


1 «2_ 


J(C,ejy+(^-u)vf 


*ta.n  — 

*  u)f  -  lO7- 


(1.27) 


Substitution  of  (1.26)  into  (1.15)  and  (1.16)  yields 

c 6  +  A,  cos>(uyfc-c*  =  0 

r  <pV'  i  *  fc  2>J  (1.28) 

8  +CqB  +•  [u£  -  M  ccs  (u)t  -c* ^  0  -  f%0si*\  tot 

which,  in  terms  of  a  new  time  variable  "C  defined  by  'll  =  6jt  -  c*.  ^  , 
may  be  written  as 

dV 


jt»+2'S.&*(Vvso*t^*0 

t—  +2  +  (S*  -t-fca  cos  O  0  =  sw  «.+<o 


where 


2#C* 


dt1 


OJ  5 


(1.29) 


*dtr 


'&TCe 


c  -  <4 
-  “73^  > 


* 


=  K*j  (1.30) 


oO" 


2V^ 


S6 

6  cO 


cO1 


(1.31) 


Both  equations  (1.29)  are  in  standard  form  except  for  the  inhomogeneous 
term  on  the  right  hand  side  of  the  equation  for  0  However,  as 
Hsu  points  out  in  [3],  for  such  an  equation  in  which  the  inhomogeneous 
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term  is  harmonic  with  a  frequency  equal  to  that  of  the  periodic  co¬ 
efficient  appearing  in  the  equation,  the  necessary  and  sufficient 
condition  that  solutions  be  stable  (unstable)  is  that  the  solutions 
to  the  corresponding  homogeneous  equation  be  stable  (unstable) .  For 
a  detailed  discussion  of  this  see  [8]. 

We  are  therefore  led  to 

Result  1:  Under  the  assumptions  leading  to  equations  (1.14),  (1.15), 
and  (1.16),  the  stability  of  the  pitching  and  rolling 
motions  of  a  ship  moving  forward  with  uniform  velocity 
in  longitudinal  waves  is  determined  by  the  positions  of 
the  pair  of  points  <  ^  •  S  >  and  (  £>$  >  )  in  the 

Ince-Strutt  chart,  where  %  .  ,  6^,  ,  $ q  ,  and  are 
given  by  (1.20)  and  (1.30).  Instabilities  are  most 
likely  to  occur  when  is  near  1/4. 

A  detailed  analysis  of  the  stability  of  the  solutions  of 
equations  (1.19)  is  too  lengthy  to  be  included  here.  Hence  we  shall 
only  sketch  the  method  used  by  Hsu  [4]  and  state  the  result. 

By  denoting  the  solutions  of  (1.19)  by  ^(.b)  and  6* (b)  and  by 
introducing  perturbations  of  these  solutions,  denoted  by  X,  and 
respectively,  one  arrives  at  the  equations  of  the  first  variation 

that 

—  :U}^  and  that  ^bt)  and  6* kb)  may  be  approximated  with  sufficient 
accuracy  by  ;*(£!)= i*  C.o%ujt  and  -  B0 Cos  LJ~t  ,  these  variational 
equations  take  the  form  of  a  set  of  two  coupled  Mathieu  equations. 
Hsu  then  shows  that  by  introducing  new  variables  ^  and  ^  related 
to  X,  and  Xj_  by  the  linear  transformation 


with  respect  to  and  6*  (b)  .  By  further  assuming 
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*  =  M  x .  -  2?  x, 

16  '  WL 


these  equations  decouple  into  two  Mathieu  equations,  the  solutions 
of  which  are  always  unstable.  The  conclusion  is 
Result  II:  Under  the  assumptions  leading  to  equation  (1.19),  the 
heaving  and  pitching  motion  executed  by  a  ship  in  calm 
water  are  always  unstable. 

We  now  turn  to  equations  (1.20).  Introducing  the  notation 

vi- 


they  can  be  written  as 

e  +.  e  »  o 


(1.32) 


The  first  of  (1.32)  can  be  Integrated  to  give  8  «d0Goe(td^t+QY  where  Q0 
and  ^  are  constants  of  integration  and,  without  loss  of  generality, 
we  may  set  Q  =  O  Substituting  the  resulting  expression  for  0 
into  the  second  of  (1.32)  gives  the  Mathieu  equation 


0eCoS  *0 


(1.33) 


which  can  be  put  into  standard  form  by  the  change  of  independent 
variable  defined  by  •£*=  60gt  .  The  result  is 

4V* 


cV-o 


d< 


where 


2/C  = 


_-/C 


oJ 


6 


oJ; 


r  _  ^4 a  A 

* 


(1.34) 


We  ere  therefore  led  to 


Eeault  HI:  Under  the  assunptions  leading  to  equation  (1.20)  the 
stability  of  the  rolling  motion  of  a  ship  oscillating 
in  pitch  and  roll  in  calm  water  is  determined  by  the 
position  of  the  point  (  ^  )  in  the  Ince-Strutt 

chart  where  ^  and  6^  are  given  by  (1.34).  Insta¬ 
bilities  are  most  likely  to  occur  when  is  near  1/4. 

Equation  (1.23)  has,  in  effect,  already  been  discussed  since  it 
is  identical  with  (1.19),  and  the  discussion  of  (1.21)  and  (1.22) 
proceeds  in  a  manner  similar  to  that  for  (1.20),  the  only  difference 
being  the  absence  of  damping  terms  in  the  former.  This  then  com¬ 
pletes  the  discussion  and  we  return  to  equations  (1.10)  (1.11)  and 
(1.12),  the  Integration  of  which  is  the  main  subject  of  this  paper. 

Equations  (1.10),  (1.11),  and  (1.12)  constitute  a  set  of  three 
nonlinear,  coupled,  second  order  differential  equations,  and  con¬ 
sidering  the  present  status  of  knowledge  concerning  such  equations, 
it  la  impossible  to  find  their  general  solutions  by  analytical  means. 
Therefore,  particular  solutions  will  be  obtained  by  means  of  an 
electronic  analogue  computer  under  the  conditions  that  the  ship  is 
initially  at  rest.  However,  before  proceeding  to  that  phase  of  the 
study  it  will  be  advantageous  to  gain  some  insight  into  the  problem 
by  making  some  simplifying  assumptions  in  the  equations  of  motion 
and  proceeding  analytically  as  far  as  possible,  just  as  was  done  in 
all  preceadlng  analyses.  We  will  show  that  under  such  assumptions 
the  equations  of  motion  for  pitch  and  roll  reduce  to  the  form  of 
inhomogeneous  Mathieu  equations,  the  stability  theory  of  which  is 
well  established.  It  should  be  stressed,  however,  that  is  is  not 
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our  Intention  to  test  the  validity  of  these  assumptions  by  means  of 
the  analogue  simulation  which  is  described  later.  Instead,  the 
analysis  which  follows  is  only  intended  to  aid  the  analogue  compu¬ 
tation  by  Indicating  if  and  under  what  conditions  instabilities  may 
be  expected. 


2.  Simplified  Analysis  Based  on  Mathieu's  Equation. 


A  great  simplification  in  the  equations  of  motion  occurs  when 
the  ship  la  assumed  to  be  symmetrical  fore  and  aft.  Under  this 
assumption 


i  -M,  =2  =  2.  =  2..  =  o 
^  =  =M;  =M--  =0 

and  equations  (1.10)  through  (1.12)  become 


(2.1) 

6  »  6  +  0  +  (t) 

(2.2) 

&  9 +  i  + 

(2.3) 

Moreover,  we  neglect  the  tens  i  ^  ^  end^|^8l 

in  (2.1) 

for  the  following  reason i  luppoee  chat  the  pitching  end  rolling 
aotlone  eve  known  functions  of  time  end  It  ie  required  to  find  ^ 
lubatituting  these  known  functions  for  ^  end  6  into  (2.1)  would 
■eke  the  tens  end  known  funotions  of  tiae 

which  could  then  be  combined  with  3(4)  ,  the  result  being,  es  fer 
eo  ^  is  concerned,  e  modification  in  the  externel  excitation,  end 
it  is  not  likely  thet  such  e  modification  would  greetly  eltev  the 
phenooenon.  Mote,  however,  thet  e  aimilar  ergunent  oould  not  be 
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used  to  neglect  the  terms  and  in  (2-2>  *nd  (2-3). 

For  suppose  ^  is  a  known  function  of  time  and  It  Is  required  to 
solve  for  p 5  and  &  Substituting  this  known  function  for  ^  Into 
(2.2)  and  (2.3)  would  have  the  effect  of  Introducing  time  dependent 
coefficients  Into  the  equations,  which  In  turn  would  greatly  alter 
the  character  of  the  solutions. 

For  this  analysis  we  will  also  disregard  the  nonlinear  damping 
factor  In  roll.  With  these  simplifications  then  (2.1)  through  (2.3) 
become 


(2.4) 

(2.5) 

*♦<«$-*>  jlsV=<^ 

(2.6) 

where  we  have  put 

=  = 

0)^  —  2,  ,  S* 

to  conform  with  the  usual  notation  used  In  the  theory  of  vibrations. 

Assusing  sinusoidal  waves,  the  external  excitation 

terms  , 

M(t")  ,  end  K("t)  may  be  written  in  the  form 

z?60* 2,  CoS  co"t  4  22SUv  Ld"t 

(2.7) 

M (t) *  M ,  cos  cot  4  Mz  siv\ 

(2.8) 

<<,  CoS  cot  4  k2siv\  U)t 

(2.9) 

or 

£(t)  =  Z?0  COS  (.  U)  "t  4  $ •s') 

(2.10) 
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M(0*  M0  ccfe  Cu)t  +  td) 


(2.11) 

(2.12) 


KCt)  *  K0cos  Cojt+tf) 
where  u)  is  the  circular  frequency  of  wave  encounter  which  depends 
on  the  geometry  of  the  waves,  the  heading  of  the  ship  with  respect 
to  the  direction  of  wave  travel,  and  the  speed  of  the  ship.  When 
(2.10)  Is  substituted  into  (2.4)  it  can  be  solved  for  7^  ,  giving 


A^cvb  (.uyfc+oi^') 


(2.13) 


where 

A  2 


Now  (2.13)  together  with  (2.11)  and  (2.12)  may  be  substituted  into 
(2.3)  and  (2.6)  yielding  the  Inhomogeneous  Mathiau  equations 

ti+CBe+  [cJq-  M^A^cos(wtfot^]e  ®  M0Cos(i <2.14) 

£  *  £ +  t60* - (cot -t-d^Sl 0  =  /Co c°^> (uit f ^ )  (2.15) 

These  can  be  put  into  standard  form  by  the  change  in  time  scale  de¬ 
fined  by  “C  *60*t’  +  o(^  .  The  result  is 

+ 2 ^  jx  ^ 6  = 

J*.  +Zk*d£  +((*t  +^cas‘C)?*feo*t‘C-dh+fy) 


(2.16) 

(2.17) 


where 


(2.18) 


*«  =  ztr  >  “a  -  ^  > 
r/=tS:  >  »  S  ”  - 

Inhomogeneous  Mathleu  equations  such  as  these  where  the  Inhomo 
geneous  term  has  the  same  frequency  as  the  periodic  coefficient 
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appearing  In  the  equation  have  already  been  discussed  above.  There 
it  was  stated  that  the  stability  of  solutions  of  such  equations  is 
the  same  as  the  stability  of  the  solutions  of  the  corresponding 
homogeneous  equations,  the  stability  chart  of  which  is  given  in  Fig¬ 
ures  1  and  2.  From  it  we  conclude  that  for  large  values  of  damping, 
unstable  solutions  are  not  likely  to  occur.  On  the  other  hand,  for 
small  values  of  damping,  unstable  solutions  are  most  likely  to  occur 
when  ^  is  in  the  neighborhood  of  1/4  and  1.  Since  damping  in  pitch 
is  of  the  order  of  50%  of  critical,  it  is  not  likely  that  unstable 
pitching  motion  will  occur.  However,  since  damping  in  roll  is  only 
of  the  order  of  57.  of  critical  (assuming  no  artificial  devices  such 
as  bilge  keels  are  used) ,  based  on  the  abov^  remarks  we  might  expect 
large  roll  angles  for  values  of  cO  near  2  03^  and 

With  this  background  we  are  now  in  a  better  position  to  investi¬ 
gate  the  solutions  of  (1.10)  through  (1.12)  on  the  analogue  computer. 
However,  before  we  can  perform  the  simulation  we  must  choose  definite 
values  for  the  coefficients  appearing  in  the  equations,  which  is  the 
subject  of  the  next  section. 

3.  Determination  of  Coefficients. 

Since  it  would  vs  prohibitive  to  examine  solutions  of  (1.10) 
through  (1.12)  for  a  variety  of  ship  forms,  we  shall,  for  concrete¬ 
ness,  adopt  the  parent  form  of  the  Series  60,  CB  -  0.60,  5'  ship 
model.  Many  of  the  coefficients  appearing  in  these  equations  may 
be  found  from  the  lines  drawing  of  the  model,  the  values  of  which 
are  given  on  the  next  page. 
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The  calculations  leading  to  these  values  may  be  found  in  the  Appendix. 
The  remaining  coefficients  are  known  to  depend  on  the  speed  and  fre¬ 
quency  of  oscillation  of  the  model,  and  curves  for  determining  these 
terms  may  be  f ound  in  [9] .  However ,  for  mathematical  convenience 
we  will  take  these  coefficients  as  constants  equal  to  a  value  within 
their  range  of  variation.  In  this  manner  we  find 

kvi-I..  *1.85  lb*4tcV-f+.  I  =  2.Z  Ib-U-sec"* 

2.  =-*>.5-  Ifc-fb-sa-c 

£.  =  4.3  lb-sec.  lb-sec. 

e  d 

£..  =  0  FL=o 

6  l 

Curves  were  not  readily  available  for  the  determination  of  the 
apparent  mass  moment  of  inertia  and  damping  coefficients  appearing 
in  the  rolling  equation  of  motion.  Concerning  the  former,  however, 
a  value  has  been  found  experimentally  for  this  model  (See  [5])  for 
rolling  in  calm  water.  It  is 

I  -  io  =  0.0454  sec7" 

*  ^ 
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and  will  be  ueed  here.  Concerning  the  letter,  we  will  take  values 
for  the  damping  coefficients  equal  to  those  obtained  from  the  "curve 
of  extinction"  as  outlined  In  [10]  using  a  free  oscillation  record 
of  the  rolling  motion  In  calm  water  obtained  experimentally  (See  [5]). 
We  thus  find 


k^  =  -O.OIS8 

Ib-^-sec  .  -0.0342  llo-^-Sec.’" 

* 

Using  these  values, 

we  then  find  from  (1.6) 

2^*  sec* 

=  *86.0  sec2- 

K^« -2.4.4  sec7" 

H;  ■*  -5. if!  s*c 
o 

-  -2.5” *1 

K*s  ■»  -  o. 348  sec- 1 

E4  =  -7.<JS'  */ie- 

—  -4>.fc8  ^-scc** 

M^=-0.5’£70  <ff-sec. 

K  =-244  Ht-stc1- 

-1.83  sac7' 

c  se 

~M  --5-4.6  &'-secv 

0  0 

2  *  -130  sec2 

M  =  -288  ^+-s«c: 

1 fi 

The  only  terms  remaining  unspecified  are  the  terms  representing 
the  external  excitation.  For  a  given  model  they  will  depend.  In 
general,  on  the  geometry  of  the  waves  and  the  heading  of  the  model 
with  respect  to  the  direction  of  wave  travel,  denoted  here  by  oC  . 
We  shall  measure  o<  In  such  a  way  that  ol»0  corresponds  to  head 
seas,  o( ■  90°  corresponds  to  waves  coming  from  the  starboard  beam, 
and  o<  -  180°  corresponds  to  following  seas. 
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From  the  •lmplified  analysis  baaed  on  Mathleu's  equation  we  may 
expect  that  for  a  given  wave  length,  X  ,  and  model  heading,  ot  , 
there  will  be  certain  values  of  u)  ,  or  equivalently,  certain  values 
of  the  models  speed,  UL  ,  for  which  the  rolling  motion  will  be  un¬ 
stable.  We  shall  seek  to  determine  this  speed  range  for  the  following 
three  cases:  (i)  X  is  taken  equal  to  L,  the  model  length,  and 

OC  assumes  the  values  15°,  30°,  45°,  60°,  75°,  105°,  120°,  135°, 
150°,  165°;  (ii)  X  is  taken  equal  to  1.5  L  and  c*  assumes  the 
values  15°,  45°,  75°,  105°,  135°,  165°;  (ill)  X  is  taken  equal  to 
0.75  L  and  0(  assumes  the  same  values  as  in  case  (11).  In  all 
cases  the  wave  height,  2A,  is  taken  equal  to  ^/40. 

The  external  heaving  force  and  rolling  and  pitching  moments 
produced  by  the  waves  were  computed  using  the  Froude -Krylov  theory 
in  conjunction  with  the  "long  wave  approximation."  The  formulas  used 
for  these  computations  are  presented  in  the  Appendix  and  the  values 
obtained  are  given  in  Table  I.  However,  before  these  values  may  be 
substituted  in  (2.7)  through  (2.9)  they  must  be  divided  by  the 
apperent  mass  of  the  model  in  the  case  of  heaving  and  the  corre¬ 
sponding  apparent  mass  moments  of  inertia  in  the  case  of  pitching 
and  rolling. 

Using  the  above  values  for  the  coefficients  the  equations  of 
motion  (1.10)  through  (1.12)  become 


^~-7q,4  y  vsi  ^-5.0,da-ni  61--  130  ^8 

-745“  6  + 

g  =  -88.0  6 -2.5^  8-288  ^6 -1.6896*  SH.  8 
-8.88  -yo.540 

£--Z4.4*d  -0.346£-0,754£\#|“244$»^ 

- 17 7  B  +  < ,  cit  Lot 


(3.1) 


(3.2) 


(3.3) 
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Table  I. 

Case  (1):  A 

-  L  ,  2A  -  ^740 

o'  -  15°: 

<i  - 

0.0172  lb-ft. 

M,  - 

4.95  lb-ft. 

f,  -  0.610  lb. 

<2- 

0.0122  " 

M2- 

0.437  " 

Z1-  -1.99  " 

ot  -  30°: 

R.  - 

0.0350  lb-ft. 

M.- 

5.25  lb-ft. 

2,  -  0.590  lb. 

<1- 

0.0230  " 

m2- 

0.333  " 

2^-  “2.55  " 

o<  -  45° 

R.  - 

0.0535  lb-ft. 

M,- 

5.25  lb-ft. 

2,  -  0.615  lb. 

*1" 

0.0324  " 

Rl- 

0.056  " 

-3.64  " 

ot  -  60° 

K,- 

0.0685  lb-ft. 

Mi  - 

4.48  lb-ft. 

2,  -  0.555  lb. 

*1- 

0.0326  " 

hr 

-0.460  " 

2L-  -  4.95  " 

o(  -  75° 

0.0825  lb-ft. 

M,- 

3.02  lb-ft. 

2,-  0.374  lb. 

Rj.- 

0.0243  " 

-0.830  " 

2^  -5.95  " 

O*  -105° 

IC,  - 

0.0825  lb-ft. 

M,- 

-3.02  lb-ft. 

2,  -  -0.374  lb 

Ki- 

-0.0243  " 

Mi" 

-0.830  " 

Z**  "5.95  " 

-120° 

R,  - 

0.0685  lb-ft. 

R,- 

-4.48  lb-ft. 

2,  1  -0.555  lb 

^2." 

-0.0326  " 

Ri- 

-0.460  " 

-4.95  " 

Ot  -135° 

K,  - 

0.0535  lb-ft. 

Mi  - 

-5.25  lb-ft. 

2,  m  -0.615  lb 

Ki." 

-0.0324  " 

K- 

0.056  " 

2^-  -3.64  ” 

O* -150° 

K.- 

0.0350  lb-ft. 

M|  - 

-5.25  lb-ft. 

2,-  -0.590  lb 

-0.0230  " 

Mi* 

0.333  " 

iir  “2-55  " 

o<  -165° 

0.0172  lb-ft. 

Mi  - 

-4.95  lb-ft. 

2,-  -0.610  lb 

-0.0122  " 

M2" 

0.437  " 

2^  -1.99  " 
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Table  I. 


Case  (1):  X  m  L  ,  2A  ■  ^/40 


ex'  -  15°: 

K,  - 

0.0172 

lb-ft. 

M,- 

4.95  lb-ft. 

2,  -  0.610 

lb. 

0.0122 

91 

Hz" 

0.437  " 

21-  -1.99 

II 

<X  -  30°: 

K.  - 

0.0350 

lb-ft. 

M,- 

5.25  lb-ft. 

2,  -  0.590 

lb. 

*z- 

0.0230 

ii 

m2- 

0.333  " 

-2.55 

fl 

o<  -  45° 

fc. - 

0.0535 

lb-ft. 

M,- 

5.25  lb-ft. 

-  0.615 

lb. 

*z- 

0.0324 

ii 

R,- 

0.056  " 

-3.64 

H 

o^-60° 

K,- 

0.0685 

lb-ft. 

4.48  lb-ft. 

2,-  0.555 

lb. 

^z" 

0.0326 

•i 

Mz" 

-0.460  " 

-  4.95 

OL  -  75° 

<•- 

0.0825 

lb-ft. 

M.- 

3.02  lb-ft. 

I,  -  0.374 

lb. 

<z." 

0.0243 

n 

Hr- 

-0.830  " 

Zj!"  -5.95 

it 

0<  -105° 

IC,  ■ 

0.0825 

lb-ft. 

M,- 

-3.02  lb-ft. 

2,  r  -0.374 

lb 

^Z- 

-0.0243 

" 

Mz- 

-0.830  " 

-5.95 

fi 

-120° 

0.0685 

lb-ft. 

M,- 

-4.48  lb-ft. 

I,  *  -0.555 

lb 

^2." 

-0.0326  " 

Mz- 

-0.460  " 

-4.95 

VI 

o(  -135° 

Kt- 

0.0535 

lb-ft. 

M,- 

-5.25  lb-ft. 

2,  «»  -0.615 

lb 

<z“ 

-0.0324  " 

PlL- 

0.056  " 

-3.64 

vi 

O* -150° 

K(  - 

0.0350 

lb-ft. 

M,- 

-5.25  lb-ft. 

2,-  -0.590 

lb 

<z- 

-0.0230  " 

Mz- 

0.333  " 

lx.-  -2.55 

ii 

o(  -165° 

K<- 

0.0172 

lb-ft. 

M,- 

-4.95  lb-ft. 

2,  -  -0.610 

lb 

Kj,» 

-0.0122  " 

Mz- 

0.437  " 

2x"  ‘1-99 

19 

1* 

1 

tfb 

**?»? 1 

(Con '  t.) 

l- 

5L,  2A  -  A /AO 

oi  -  15° 

a 

->  0,0202  lb-ft. 

M(-  8.48  ik^ft. 

f,-  0.913  lb. 

- 

a. 00121  " 

F\-  -0.236  " 

hm  ‘V5  " 

o(-  45? 

0.0580  lb-ft. 

■.  VI 

0.0254  " 

(jj,-  7.18  lb-ft. 
Mz-  -0.839  " 

2,  -  0.f$5  lb. 

2*-  -8.77  " 

o<  -  7si| 

K, -  0.0828  lb-ft. 

M, -  3.03  lb-ft. 

?i-  o?  0^4  .f  . 

|<:2-  0.0143  " 

M2-  -1.49  " 

2t-  -10.8  " 

o(  -105°; 

K\  -  0.0828  lb-ft. 

Mi-  -3.03  lb-ft. 

2,-  -0.334  lb. 

^  -0.8143  " 

i' 

Mr-  -1.49  " 

2Z-  -10.8 

o<  -135° 

1^VX);0580  lb-ft. 

M,  -  -7.18  lb-ft. 

2,  -  -0.795  lb. 

fi*-  -0.0254  " 

Pt2-  -0.839  " 

Z*.-  >8.77  " 

o(  -165° 

|C,  -  0.0202  lb-ft. 

M,  -  *8.48  lb-ft. 

2,  -  -0.913  lb. 

*0.00121  " 

Mz-  -0.236  " 

2t-  -7.0*  " 

dll):  A- 

0.75L,  2A  -  */40 

o(  -  15° 

0.0126  lb-ft. 

M,  -  2.28  lb-ft. 

2,  -  0.210  lb. 

|C>  0.T)0906  " 

Mz-  0.612  " 

-0.0840  lb. 

Of  -  45° 

lc,-*a.o^i  ib-ft. 

M(  -  3.36  lb-ft. 

2,  -  0.411  lb. 

0.0322  " 

F^-  0.339  " 

2Z-  -1.26  " 

o(  -  75° 

fc,-  0.0861  lb-ft. 

R,-  2.28  lb-ft. 

0.308  lb. 

0.0260  " 

**  v* 

R2-  -0.572  " 

Z2-  -3.74  " 

\ 

6  '■%  •' 
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Table  I 

(Con  1 1 . ) 

oi  -105° 

0.0861  lb-ft. 

Pj|  -  -2.28  lb-ft. 

-  -0.308  lb. 

-0.0260  " 

-0.572  " 

2r  -3.74  " 

ok  •lift* 

K,-  0.0472  lb-ft. 

M,-  -3.36  lb-ft. 

2> -  -0.411  lb. 

-0.0322  " 

M^-  0.339  " 

-1.26  " 

o<  -165° 

Kt-  0.0126  lb-ft. 

R, -  -2.28  lb-ft. 

2,-  -0.210  lb. 

j 

-0.00906  " 

Mt-  0.612  " 

1^-  -0.0840  " 
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where  use  ha»  also  been  made  of  (2.7)  through  (2.9). 

Be fora  proceeding  to  the  explanation  of  the  analogue  simulation 
let  ua  restate  the  assumptions  which  have  been  made  In  arriving  at 
values  for  the  various  coefficients  In  the  equations  of  motion. 

1.  The  coefficients  of  the  terms  representing  hydrodynamic 
forces  and  moments  are  constants  equal  to  representative  values 
obtained  experimentally. 

2.  The  waves  are  sinusoidal. 

3.  The  magnitudes  of  the  external  forces  and  moments  are  ob¬ 
tained  from  the  Froude-Krylov  theory  In  conjunction  with  the  "long 
wave  approximation." 

All  of  these  assuaptlons  seem  justified  since  It  is  not  our 
intention  to  give  a  detailed  study  of  the  motion  of  the  Series  60, 
Cg  -  0.60  ship  model.  Bather  we  are  concerned  with  studying  a 
basic  phenomenon,  and  this  model  has  been  adopted  only  to  errlve  at 
realistic,  if  not  completely  accurate,  values  for  the  various  terms 
appearing  In  the  equations  of  motion. 

Mow  that  definite  values  have  been  chosen  for  the  coefficients 
in  the  equations  they  may  be  simulated  on  the  analogue  computer, 
which  is  the  subject  of  the  next  section. 

4.  Analogue  Simulation.. 

In  simulating  a  problem  such  as  this  on  the  analogue  computer 
it  is  necessary  that  the  output  of  any  operational  amplifier  not 
exceed  t  100  volts.  To  satisfy  this  requirement,  it  is  necessary 
to  4-ntroduce  scale  factors  between  the  physical  variables  and  the 
machine  variables  (i.e.  voltages) based  on  the  estimated  slse  of  the 
former.  From  physical  considerations  we  may  state  that  the  largest 
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valua  of  |^|  to  be  expected  Is  1/6  ft.,  the  largest  value  of  |Q( 
to  be  expected  la  1/5  radian,  and  the  largest  value  of  |^5|  to  be 
expected  is  2  radians.  Thus  If  we  define 

B  *  ^oo  Y»5’oo0  X—So^  (4.1) 

the  machine  variables  X,  Y,  and  2  will  all  be  less  than  100  volts 
In  absolute  value.  Substituting  (4.1)  In  (3.1)  through  (3.3)  gives 

B  =  -79,4*  -5.51  2  —  onzo  X^-o.41/  Y^-o.ifeoYB 
-9,5-4  Y 4.  2.79  Y-t-fcoo^coeoJ"*:  +  k22siir>aJ"t 

V  -  -8fc,o  Y  -2.59  Y  -0.48  Y2  -o.54.fc  Xx-o.o75S  2* 

-5. 57  2  -0.492  2  4-5*00  HjCofcuit^Soo^^w  oit 

X  *  -24.4  X -0.848  X-0.OI49  X  jx| -o.4o7  X2 
-0.554  X  Y  +  5o  K^Coscjt+iTOiq^tjtr 


(4.2) 


(4.3) 


(4.4) 


It  Is  also  convenient  to  have  the  absolute  values  of  all  the  co¬ 
efficients  In  the  equations  less  than  1  which  can  be  accomplished 
by  a  change  of  time  scale  defined  by 

T-lot  (4.5) 

where  T  is  now  the  machine  time  and  i:  Is  real  time.  After  this 

change  of  variable  we  obtain  the  final  equations 


2  «  -0.794  i  -0.35-1  2 


/  0.720  v2- 


.1.  n.TLD 


-0. 09  54  Y  4  0. 2 7  9  Y  V  fc 2 ,  COS  ^  T+ tZj Si  7§  7" 
YM-0.8fc0Y-0.254Y'-^  2* 

-0.05S"7  £--0,0492  2Vs’M,GoS^'T«’5‘Hl4in^T 


(4.6) 


(4.7) 
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A/'=-0.24^X-O.OS4aX,-O.Ol4.«»X'|X'|-^2  XE 

_aj|aXy+0.SKiW^T+0-5.Ki<>m^T  (4-8) 

where  primes  denote  differentiation  with  respect  to  T. 

The  simulation  of  these  equations  is  done  in  accordance  with 
the  block  diagram  shown  in  Figure  3.  The  triangles  represent  phase 

inverters,  the  triangles  with  rectangles  along  one  side  represent 

/  4 

summing  integrators,  the  circles  on  the  input  sides  of  the  Integrators 
represent  potentiometers,  the  initial  conditions  (if  any)  are  shown 
above  the  Integrators  concerned,  and  the  multiplication  by  a  con¬ 
stant  performed  through  each  Integrator  or  phase  inverter  is  shown 
opposite  the  corresponding  input.  The  computer  used  was  a  Berkeley 
Ease  (Electronic  Analogue  Simulating  Equipment)  electronic  analogue 
computer  manufactured  by  Berkeley  Division  of  Beckman  Instruments, 

Inc.,  and  is  shown  in  Figure  4. 

The  simulation  of  a  complex  problem  such  as  this  on  an  analogue 
computer  presents  many  opportunities  for  making  errors,  particularly 
in  making  the  external  connections  into  the  patch  bay  of  the  computer, 
and  since  the  solutions  of  equations  (3.1)  through  (3.3)  are  not 
known,  such  errors  may  not  be  easily  detectable.  However,  these 
equations  contain  equations  (2.4)  through  (2.6)  as  a  special  case 
which  may  be  recovered  by  setting  certain  nonlinear  terms  equal  to 
zero.  Since  the  stability  of  these  latter  equations  is  well  estab¬ 
lished  we  may  use  that  theory  to  check  the  validity  of  the  simulation. 

To  do  this  the  patch  board  was  first  wired  to  simulate  equations 
(3.1)  through  (3.3)  in  accordance  with  the  block  diagram  shown  in 
Figure  3.  Then  appropriate  coefficients  were  set  equal  to  zero  so 
that  the  resulting  simulation  was  that  of  equations  (2.4)  through 
(2.6)  except  now  definite  values  have  been  assumed  for  the  coefficients. 
We  then  sought  the  values  of  cO  for  which  the  rolling  was  unstable. 


FIG.  3  BLOCK  DIAGRAM  FOR  EQUATIONS  (4.6),(4.7) ,  (4.8) 


FIG.  3  (CONTINUED) 


Figure  4  Berkeley  EASE  Electronic  Analogue  Computer 


-34- 


The  case  A-  L  ,  o< -  45°  was  chosen  for  the  computation  and 
samples  cf  the  rolling  response  obtained  for  values  of  60  near  the 
stability  threshold  are  shown  in  Figure  5.  It  can  be  seen  that  near 
the  threshold  a  small  change  in  OJ  produced  a  large  change  in  the 
character  of  the  solution.  The  unstable  solutions  showed  a  sus¬ 
tained  growth  in  amplitude  which  would  have  continued  indefinitely. 

It  was  found  that  the  rolling  motion  was  unstable  for  8.28<u)< 
10.5  whereas  the  theory  predicts  unstable  rolling  for  8.26<(*) <10.6. 
The  error  in  the  bounds  of  the  range  of  U)  for  which  the  rolling 
was  unstable  obtained  by  the  computer  is  less  than  1%  which  seems  to 
indicate  that  the  basic  network  was  correctly  constructed.  There¬ 
fore,  it  seems  reasonable  to  conclude  that  the  network  will  give 
the  correct  particular  integrals  of  (3.1)  through  (3.3).  The  pre¬ 
sentation  and  discussion  of  these  integrals  is  the  subject  of  the 
next  section. 

5.  Computer  Results. 

i 

It  was  found  that  for  the  cases  where  A-  L  and  A-  0.75L 
the  rolling  response  had  the  general  appearance  shown  in  Figure  6. 

It  was  found  that  for  these  two  cases  there  was  a  range  of  exciting 
frequencies  for  each  value  of  o(  for  which  the  roll  amplitude  tended 
to  grow,  just  as  was  observed  for  the  solutions  of  the  Mathleu 
equation  shown  in  Figure  5.  However,  because  of  the  second  order 
damping  term  in  the  rolling  equation  the  growth  in  amplitude  did  not 
continue  indefinitely,  but  instead  reached  a  steady-state  value. 
Nevertheless,  because  of  the  similarity  to  the  solutions  of  the 
Mathleu  equation  we  will  term  such  motion  "unstable."  To  give  some 
idea  of  the  rolling  amplitudes  attained  we  refer  to  Figure  7  which 
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Figure  b  Rolling  Response  as  Governed  by  Eqs.  (2.4  ,  (2.5'' 
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Figure  5  (cont.)  Rolling  Response  as  Governed  by  Eqs.  (  2  .  (2.5) 

A=L,  o<  =45° ,  2A=  V:0 
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(d)  Unstable  Solution, 
uJ«7.1o  rad/sec 


Figure  b  (cont.)  Rolling  Response  for  A-L,  o<«45°,  2A**  A/40 


cu,  rad/sec. 


FROUDE  No. 

FIG.?  ROLLING  AMPLITUDE  AS  A  FUNCTION 
OF  w  FOR  X'L,  a  *45*  2A*  X/40 

go,  rad /sec. 


0  0.1  0.2  0.3 


FROUDE  No. 

FIG.  8  ROLLING  AMPLITUDE  COMPUTED  FROM 
LINEAR  THEORY  FOR  \=L,  a *45®  2A= 
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shows  the  maxinun  amplitude  as  a  function  of  UJ  and  Froude  number 
for  o<«  45°  and  A*  L,  and  for  comparison  the  same  curve  is 
plotted  in  Figure  8  in  which  the  rolling  amplitude  was  calculated 
using  linear  theory. 

Although  the  rolling  response  did  exhibit  the  instabilities 
just  described  for  A-  L  and  X  -  0.75L,  .they  did  not  occur  at 
values  of  60  predicted  by  the  theory  based  on  the  inhomogeneous 
Mathieu  equation,  which  is  not  surprising  considering  the  assumptions 
made  in  that  development.  Instead  the  unstable  regions  were  found 
to  be  much  broader,  in  general,  than  those  predicted  by  that  theory, 
as  may  be  seen  by  Table  2,  which  shows  the  comparison  for  the 
special  case  A  -  L,  2A  -  A/ 40.  Nevertheless,  as  may  be  seen, 
the  instabilities  did  occur  for  values  of  in  the  neighborhood 

of  1/2  as  the  theory  predicts. 

For  the  case  where  A  -  1.5  L  the  rolling  response  was  found  to 
be  quite  different  from  that  obtained  for  A  -  L  and  X  ■  0.75  L, 
as  may  be  seen  by  Figures  9  and  10.  Because  of  the  very  large  roll 
angles  obtained  the  motions  found  in  this  case  are  very  difficult  to 
interpret  physically  since  the  equations  of  motion  are  no  longer 
valid.  Nevertheless,  the  results  are  interesting  since  they  clearly 
show  the  effect  of  the  nonlinear  terms  when  the  displacements  are 
large . 

The  solutions  obtained  showed,  in  general,  that  for  each 
heading  of  the  model  with  respect  to  the  direction  of  wave  travel 
there  were  values  of  CO  for  which  the  rolling  was  unstable.  The 
model's  speed  corresponding  to  these  values  of  was  then  found 
from  the  following  equation 


(I):  Rang*  of  for  which  solutions  to  the  Hathieu 

equation  (2.6)  are  unstable. 

(II):  Range  of  for  which  solutions  to  equation  (3.3) 

are  unstable. 

In  both  cases  A  -  L,  2A  -  ^/40. 


ire  9  Rolling  ResDonse  for  A-1.5L,  o<=l5°,  2A=  A/40 


43 


-1.0  * 
-0.5  J 


(b)  Unstable  Solution, 
uJ-3.88  rad/sec 
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(c)  Stable  Solution,  cJ-4.00  rad/sec. 


Figure  10  Rolling  Response  for  A=*1.5L,  o(=105°,  2A=  A/70 


u-=  — * —  (  rbr<^'^uj) 

Cos  oC  '  2TT  w  ' 


where  Is  the  velocity  of  gravity  waves  given  by 


The  results  of  the  stability  study,  may  be  found  In  Figures  11 
and  12  which  show,  for  each  value  of  ot  ,  the  range  of  Froude  num¬ 
bers  for  which*  the  rolling  was  found  to  be  unstable.  For  "Froude 
manber8  between  0  and  0.35  it  was  found  that  the  rolling  motion  was 
always  stable  for  values  of  o(  >  105°  for  the  cases  where  X "L  and 
A  *0.75  L,  and  hence  these  parts  of  Figure  11  are  omitted. 

In  addition  to  the  growth  in  amplitude,  the  unstable  rolling 
motions  had  the  unusual  characteristic  that  the  period  of  the  steady- 
state  response  had  a  value  equal  to  twice  the  period  of  wave  encounter, 
while  the  stable  motions  had  a  period  equal  to  the  period  of  the 
external  excitation.  This  is  clearly  seen  in  Figure  13  (a)  and  (b) . 

In  the  former  the  period  of  wave  encounter,  T  ,  corresponding  to 
t *)m  8.84  rad/sec.  is  0.711  sec.  while  from. the  motion  record  the 
period  of  roll  is  found  to  be  1.4  secs.  In  the  latter,  the  value 
of  T  corresponding  to  U)  ■  7.01  rad/sec.  is  0.897  secs,  while  from 
the  motion  record  the  period  is  found  to  be  0.9  sec.  Figure  6  (b) 
shows  an  example  of  a  rolling  response  near  the  threshold  of  the  un¬ 
stable  region  where  the  time  between  alternate  large  rolls  is  again 
equal  to  twice  the  period  of  wave  encounter.  We  will  have  more  to 
say  about  this  phenomenon  in  the  next  section. 
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FIG.  II  RANGE  OF  FROUDE  Nos  FOR  WHICH 
ROLL  IS  UNSTABLE 
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FIG.  12  RANGE  OF  FROUDE  Nos  FOR  WHICH  ROLL  IS  UNSTABLE 

X  *  I.5L  ,  2A»  V40 


6.  Discussion  and  Conclusions. 


We  have  seen  that  when  the  frequency  Of  wave  encounter,  tO  ,  is 
in  the  neighborhood  of  twice  the  natural  frequency  in  roll,  , 

the  rolling  response  has  the  following  two  distinctive  characteristics 

(I)  the  amplitude  tends  to  grow  to  a  large  but  finite  value,  and  (2) 
the  period  of  the  steady  -.state  response  is  equal  to  twi.ee  the  period 
of  the  excitation.  This  behavior  cannot  emerge  from,  nor  ean  it  be 
studied  in,  a  linear  analysis.  For  if  all  the  nonlinear  terms  are 
neglected,  the  equation  of  motion  for  roll  is  no  longer  coupled  with 
those  for  heave  and  pitch,  and  is  linear.  One  is  then  led  to  the 
conclusions  that  large  roll  amplitudes  may  be  expected  only  when  the 
exciting  frequency  is  near  the  natural  rolling  frequency, . the  reson¬ 
ance  condition,  and  that  in  any  case  the  frequency  of  the  steady- 
state  response  is  equal  to  the  frequency  of  the  excitation.  There¬ 
fore,  we  see  that  the  inclusion  of  the  second  order  static  coupling 
terms,  which  are  known  to  exist,  in  the  equations  of  motion  produces 
a  rolling  behavior  which  is  in  direct  contradiction  to  that  predicted 
by  linear  theory.  However,  in  this  study  we  af.e  dealing  only  with 
solutions  to  certain  differential  equations,  and  whether  or  not  the 

i  .  i 

rolling  motion  of  an  actual  ship  moving  in  oblique  waves  would 
exhibit  such  behavior  depends  upon  how  well  those  equations  approxi¬ 
mate  the  true  equations  of  motion,  a  question  which  can  only  be 
answered  by  model  experiments. 

Some  experiments  with  Series  60  models  moving  in  oblique  waves 
have  recently  been  performed,  the  results  of  which  may  be  found  in 

[II]  and  [12].  The  results  show  that  in  bow  seas  the  roll  amplitude 
was  always  of  the  order  of  5°  for  Froude  numbers  between  0  and  0.4, 


«  it  ..  k 
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«fhi&  doe*  not  conform  with  our  results  shown  in  Figure  7.  However, 
thediScfepancy  may  be  attributable  to  two  main  factors.  First  of 
all,  in  deriving  the  equations  of  motion  It  waa  assumed  that  all 


all,  in  deriving  the  equations  of  motion  It  waa  assumed  that  all 

^  -j.'. 

motions  are  small  in  the  sense  that  third  order  terms  may  be 
neglected,  and  therefore,  for  such  large  roll  angles  4h*  validity 
of  those  equations  is  questionable.  It  seems  apparent  freer  the  com¬ 
puter  study  that  within  a  certain  range  of  Fromde  numbers  the'rolllng 
motion  may  tend  to  become  large,  but  how  large  is  difficult  toP^say 
at  this  point.  Secondly,  the  models  used  in  the  experiments  described 

»y**  *  ^  1  . 

in  [11]  were  all  equipped  with  bilge  keels,  and) in  both  ftudles  the 

A'M,  •  "  •  u 

waves  used  in  the  experiments  were  not  very  severe.  (In  [II]  and  [12] 
the  wave  heights  were  taken  equal  to  A/ 50  and  A/48  respectively). 
Therefore,  it  seems  likely  that  the  damping  in  roll  wad  too  large 
compared  to  the  wave  excitation  to  produce  unstable  rolling. 

On  the  other  hand,  rolling  motions  in  which  the  response  had  a 
period  equal  to  twice  the  period  of  wave  eneounter  >■  were  obaerved 
during  the  experiments  reported  in  [12],  where  a  recorded  rolling 
motion  is  presented  (Figure  24  of  that  paper)  that  very  closely  re¬ 
sembles  that  shown  in  Figure  6  (b)  of  this  paper.  In  connection 
with  this,  Lewis  and  Numata  state  that  such  rolling  motions  wars 
observed  when  the  period  of  wave  encounter  was  approximately  one 

v.*r 

half  the  natural  period  in  roll.  However,  they  also  point  out  that 

*i 

in  other  cases  where  the  two  periods  did  have  this  relationship, 
normal  records  were  obtained.  No  definite  explanation  for  this  be- 
havior  is  given  there,  but  based  on  the  Computer  study  we  are  in  a 
position  to  offer  an  explanation. 

First  of  all  it  is  evident  from  this  investigation  that  the 
presence  of  the  nonlinear  static  coupling  between  heave,  pitch,  end 
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roll  Is  sufficient  (provided  the  damping  in  roll  Is  not  too  large) 
to  cause  such  behavior  and  that  no  motions  other  than  these  need  to 
be  Involved.  Secondly,  we  have  found  that  such  behavior  Is  a  prop¬ 
erty  of  an  unstable  or  nearly  unstable  rolling  motion  as  we  have 
defined.  Now,  for  unstable  rolling  to  occur  it  is  necessary  but 
not  sufficient,  that  the  frequency  of  wave  encounter  be  in  the  neigh¬ 
borhood  of  twice  the  natural  rolling  frequency.  Whether  or  not  un¬ 
stable  rolling  actually  occurs  depends  upon  the  relative  strength  of 
the  excitation  compared  to  the  damping  in  roll.  This  then  explains 
why  in  one  case  a  "normal"  rolling  behavior  and  in  another  case  an 
"abnormal"  rolling  behavior  can  be  observed  even  when  in  both  cases 
the  frequency  of  wave  encounter  is  in  the  neighborhood  of  twice  the 
natural  rolling  frequency.  If  the  rolling  is  unstable,  then  it 
exhibits  a  period  equal  to  twice  the  period  of  wave  encounter,  which 
as  a  consequence  is  nearly  equal  to  the  natural  period. 

To  be  completely  accurate,  the  "erratic"  rolling  motion  ob¬ 
tained  by  Lewis  and  Numata  does  not  satisfy  our  definition  of  unstable 

.  I 

rolling  since  the  amplitude  is  not  increasing  and  shows  no  tendency 
to  do  so  and  hence  the  above  explanation  is  not  strictly  applicable. 

In  fact  it  is  not  quite  periodic  so  it  is  meaningless  to  speak  of 
1 ts  period  as  being  twice  the  period  of  excitation.  However,  it 
does  closely  resemble  the  computer  solution  shown  in  Figure  6  (b) , 
and  therefore  we  may  discuss  it  with  reference  to  that  solution. 

As  may  be  seen  the  sblutlon  in  Figure  '6  (b)  is  classified  as 
stable  but  is  very  near  the  stability  threshold.  Furthermore,  the 
initial  portion  of  the  solution  is  almost  periodic,  the  time  between 
the  alternate  large  rolls  being  precisely  twice  the  period  of  the 
excitation.  However,  this  behavior  is  transient  and  gradually  dies 
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out  leaving  a  steady  atata  response  which  la  periodic,  tha  period 
balng  equal  to  the  parlod  of  the  excitation.  The  cloae  resemblance 
of  the  above  mentioned  rolling  motion  to  this  solution  seems  to  indi¬ 
cate  that  what  Lewis  and  N\snata  actually  observed  was  the  transient 
portion  of  a  rolling  motion  which  was  very  near  the  stability  thresh¬ 
old,  the  transient  being  almost  periodic  with  period  equal  to  twice 
the  period  of  wave  encounter  and  which  as  a  consequence  was  nearly 
equal  to  the  natural  period. 

Tha  investigation  described  in  this  paper  can  only  be  considered 
as  an  initial  step  in  the  investigation  of  the  effect  of  nonlinear 
static  coupling  on  the  rolling  stability  of  a  ship  moving  in  oblique 
waves.  However,  it  does  suggest  a  possible  area  of  study  for  those 
experipentors  engaged  in  the  testing  of  models  in  waves. 


-52- 


APPENDIX 


Formulas  for  the  calculation  of  the  coefficients  of  the  term! 
involving  the  displacements  in  the  equations  of  motion  are  presented 
in  [1].  Those  which  depend  only  on  the  geometry  of  the  submerged 
volute  are  —  A 


where  is  the  area  of  the  equilibrium  water  plane,  is  the 
x-coordlnate  of  the  center  of  flotation,-  ^  is  the  density  of  the 
water ,  and  dy/  is  the  slope  of  the  top  sides  at  the  water 

line  (positive  for  tumblehome).  The  remaining  eoefflcients  depend 
on  the  location  of  the  center  of  gravity  as  well  as  the  geometry  of 
tha  submerged  volume  and  are  given  by1 


K(4=-ASMt 


M4— 

jt, 

\l  *  h'u)  It A  V>\j  V 

-K 

r~ - 

Several  sign  errors  were  made  in  the  derivation  of  these  formulas 
in  {1]  which  have  now  been  corrected. 
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whare  1^  is  the  moment  of  inertia  of  the  .equilibrium  water  plane 
about  the  y-axis,  zfi  ia  the  z-coordinate  of  the  center  of  buoyancy, 
zw  is  the  z-coordinate  of  the  point  where  the  z-axis  cute  the  equil¬ 
ibrium  water  plane,  and  b(x)  ie  the  half  breadth  of  the  equilibrium 

,ir 

water  plane. 

For  these  computations  we  chose  GM^  to  be  5%  of  the  beam  or 
0.40".  From  the  geometry  of  the  submerged  volume  we  found  for  this 
model  that  BHp  ■  1.51"  and  that  the  center  of  buoyancy  was  1.46" 
from  the  water  plane.  Thus 

zw  -  1.51"  -  (0.40"  +  1.46")  -  -0.35"  -  -0.0292’ 

and 

zB  -  -0.35"  +  1.46"  -  1.11"  -  0.0925 ' 


Now  all  that  remains  to  evaluate  these  coefficients  is  to  perform 
the  numerical  Integrations. 

The  heaving  force  and  the  rolling  and  pitching  moments  pro¬ 
duced  by  the  waves  were  calculated  from  the  following  formulas 


-A^  \  '[ZbliVkSUflijcos  cix 

-i 

^  [ x \2 box'* -k «»uS]tns  kt  X  -*-k,siWkvx 

R,  ■ -  M  f  \%  i  V  yj^***, *  ix 

u  -1*. 

-  Rv  aw 


where 


S(x)  -  station  area  z,  -  vertical  offset  of 

the  base  line 

Kx  -  K  cosoC 

A  ■  wave  amplitude 

K2  -  K  sinoL 
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